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We connect the rare fluctuations of an Equilibrium (EQ) process and the typical fluctuations 
of a nonequilibrium (NE) stationary process. In the framework of large deviation theory, this 
observation allows us to introduce NE thermodynamic potentials. For continuous-time Markov 
chains, we identify the relevant pairs of conjugated variables and propose two NE ensembles: one 
with fixed dynamics and fluctuating time-averaged variables, and another with fixed time-averaged 
variables, but a fluctuating dynamics. Accordingly, we show that NE processes are equivalent to 
conditioned EQ processes ensuring that NE potentials are Legendre dual. We find a variational 
principle satisfied by the NE potentials that reach their maximum in the NE stationary state and 
whose first derivatives produce the NE equations of state, and second derivatives produce the NE 
Maxwell relations generalizing the Onsager reciprocity relations. 


I. INTRODUCTION 

Potentials define a specific concept in physics. They 
predict the evolution of a system from a variational prin¬ 
ciple. Such principles span many scientific fields from 
mechanics, electromagnetism and optics to control the¬ 
ory, thermodynamics and statistical physics. A varia¬ 
tional principle elegantly summarizes the method used to 
solve a problem into the extremization of the appropri¬ 
ated cost function, for instance the action in mechanics 
[1] , the optical path length in optics [2] , or the thermody¬ 
namic potential in statistical physics [3] . The underlying 
idea is to explore all possibilities, including non physical 
ones, to find the physical solution from the extremum of 
the cost function. 

In statistical physics, a thermodynamic potential is a 
state function of the thermodynamic variables. The lat¬ 
ter specify a coarse-grained representation of the state of 
a system including a large number of degrees of freedom. 
Thermodynamic variables come in conjugated pairs: in 
each pair, one variable is free and one is constrained ac¬ 
cording to the environmental conditions. The EQ ther¬ 
modynamic potentials proceed from the Legendre trans¬ 
formation of either energy or entropy. This transforma¬ 
tion, at the core of the theory’s dual structure, allows us 
to interchange the free and constrained variables. The 
thermodynamic state is reached at the extremum of the 
thermodynamic potential. There, the mean free variables 
are functions of the constrained ones. Beyond the mean 
description, the potential also predicts the statistics of 
the free thermodynamic variables, either by generating 
their cumulants, or from its connection with the asymp¬ 
totic probability of the free variables. 

Statistical physics provides a microscopic foundation 
to thermodynamics and a method to describe equilib¬ 
rium systems. In the last decades, the large deviation 
theory [4, 5] has modernized our understanding of sta¬ 
tistical physics and accounted for its successes. More 
recently, it has received a growing interest thanks to its 
applications to NE systems, for instance in glasses [6-9], 
biological systems [10-12] or rare events sampling [13]. 


Clearly, one step toward understanding NE phenomena 
starts with the derivation of a NE thermodynamic po¬ 
tential verifying most of the aforementioned properties. 
With this in mind, many authors have shed light on the 
structure of statistical physics for NE Markov processes. 
Oono and Eyink considered that Large Deviation Func¬ 
tions (LDF) could represent NE potentials [4, 14, 15]. On 
this basis, Oono and Paniconi proposed a phenomeno¬ 
logical framework to study NE steady states [16]. For 
NE continuous processes, Bertini et al. developed the 
macroscopic fluctuation theory describing the statistics 
of density and current fluctuations in Non-Equilibrium 
Stationary States (NESS) [17, 18]. Bodineau and Derrida 
used an additivity principle to predict those fluctuations 
in diffusive systems [19, 20]. 

For discrete processes, Lecomte, Appert-Rolland and 
van Wijland introduced a dynamical partition function 
and the corresponding topological pressure identihed as 
a LDF [21, 22]. Baule and Evans explored these ideas 
using a path entropy with the aim of hnding rules con¬ 
straining the dynamics of fluids under continuous shear 
[23-25]. Monthus proposed a similar approach, but in¬ 
volving the maximization of a trajectory-based relative 
entropy in the presence of constraints [26]. Using large 
deviation theory, Maes and Netocny [27] found a canon¬ 
ical structure and obtained the LDF of occupation and 
current probabilities from a variational approach based 
on the LDF of occupation and transition probabilities. 
A key step was the introduction of an EQ reference pro¬ 
cess to highlight that EQ fluctuations naturally appear 
when studying NE fluctuations. From another perspec¬ 
tive, Nemoto and Sasa have shown that a Cumulant Gen¬ 
erating Function (GGF) also proceeds from a variational 
principle, strengthening the dual structure of the theory 
[28]. 


More recently, Chetrite and Touchette proposed a gen¬ 
eral framework for both continuous and discrete pro¬ 
cesses: they found that a conditioned Markov process is 
ensemble-equivalent to a condition-free process called the 
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Table 1. Relationship between the various stochastic processes and NE ensembles. The EQ reference process conditioned on the 
energy currents j, activities / and occupations p generates the trajectories of the systems in the NE micro-canonical ensemble. 
The process with mean energy currents, activities and occupations equal to the constrained values of the conditioned process 
is the driven process. The path probabilities of the driven process are asymptotically equivalent to the path probability of the 
NE process and of the canonical reference process. The CGE of j, f and p for the NE process is exactly the same as the CGE 
of the EQ reference process (up to a translation), i.e spontaneous rare fiuctutions of the EQ process are associated to typical 
realizations of the NE process. The NE process generates the trajectories of the systems in the meta-canonical ensemble. This 
ensemble includes the systems that are put out-of-equilibrium by gradients of temperatures imposed by heat reservoirs. From 
the equivalence between the conditioned reference process with the driven reference process and the NE process, we conclude 
that the NE micro-canonical and meta-canonical ensembles are equivalent. 


Ensembles of NE systems 
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Similarity transformation K = \Tr\ ■ k ■ In] ^ 


Asymptotic path equivalence 


driven (or auxiliary) process, but also to an exponentially 
tilted process called the canonical process [29, 30]. This 
later process is defined by exponentially weighting the 
probability of each trajectory with a weight depending on 
a functional v of the stochastic process. This weighting 
procedure, is analogous to the definition of the canonical 
ensemble from the superposition of micro-canonical en¬ 
sembles using a Boltzmann weight. On the other hand, 
the conditioned Markov process assumes that the vari¬ 
able V is constrained to a given value. Finally, the driven 
process has a dynamics defined such that the mean value 
of V is equal to the imposed value in the conditioned pro¬ 
cess. A systematic method of constructing this driven 
process from a variational approach was provided in Ref. 
[31]. A construction of the canonical process was also 
proposed by Giardina, Kurchan and Peliti in Ref. [32] 
for classical systems and by Garrahan and Lesanovsky 
in Ref. [33] for dissipative quantum systems. Jack and 
Sollich constructed a driven process for classical systems 
in Ref. [34] . The questions of the validity of the path en¬ 
semble equivalence has recently been studied in Ref. [35] 

Despite all these results, the structure of NE statistical 
physics is incomplete as regards to EQ statistical physics. 
For instance, the identification of the relevant coarse¬ 
grained degrees of freedom, i.e. the NE thermodynamic 
variables, is still missing. Accordingly, no general defini¬ 
tion exists for stationary NE thermodynamic potentials. 
To progress in this direction, focusing on continuous-time 


Markov chains and stationary processes, we consider the 
following questions: can we describe the NE fluctuations 
of a system from the fluctuations of the same system at 
EQ? If yes, can we define meaningful NE thermodynamic 
potentials using the variables involved in this correspon¬ 
dence? We positively answer these two questions by find¬ 
ing an exact mapping between the statistics of EQ and 
NE processes. This mapping involves, among others, the 
affinities of the NE process and some dynamical biases. 
The later parameter enables to dilate the energy barriers 
separating the various states of the system. The variables 
conjugated to the affinities and the dynamical biases are, 
respectively, the energy currents and the activities of the 
exchanges with the environment. The existence of a sim¬ 
ple mapping when considering the appropriated couples 
of conjugated variables suggests that a complete canon¬ 
ical structure for NE statistical physics exists. With 
respect to previous works on conditioned Markov pro¬ 
cesses, our main contribution is to identify the constrains 
that does not modify the system dynamics, apart from 
changing the temperatures of the heat reservoirs. Ac¬ 
cordingly, we define two ensembles of NE systems: the 
meta-canonical ensemble where the constrained variables 
are the affinities, and the NE micro-canonical ensemble 
where the constrained variables are the energy currents. 
We prove the equivalence of these ensembles and derive 
the NE thermodynamic potentials conjugated by Legen¬ 
dre transformation. We also obtain the NE equations of 
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state connecting the conjugated variables. 


Our results and the structure of the theory are sum¬ 
marized in Table I. Accordingly, the outline of the paper 
is as follows. We start by studying the fluctuations of 
an EQ reference process in Sec. II whose material cor¬ 
responds to the middle row of Table I. The definition of 
the EQ reference process and an introduction to large de¬ 
viation theory are provided in Secs. IIA and IIB. After 
these introductory subsections, we look for an asymptotic 
approximation of the probability of the energy currents, 
activities and occupations of the systems states. Since 
we are dealing with an EQ system, no mean energy cur¬ 
rent exists. However, rare spontaneous fluctuations may 
produce non-zero energy currents and some arbitrary ac¬ 
tivities and occupations. We seek the probability of these 
events from an optimization problem: given that some 
energy currents j, activities /, and occupations p are ob¬ 
served, defining the conditioned reference process, which 
process (called the driven reference process) reproduces 
these conditioned values j, / and p as typical values? We 
construct this driven process in Sec. IIC and obtain the 
LDE of j, / and p. We use this result to derive the cor¬ 
responding scaled CGF from a variational approach in 
Sec. IID. 


We switch to the study of the fluctuations of a NE 
process in Sec. III. This section corresponds to the third 
row of Table I, which is obtained following exactly the 
same path as for the EQ reference process, except that 
we start with a NE process as defined in Sec. Ill A: we 
look for the NE driven process that will typically repro¬ 
duce the arbitrary energy currents j and activities / im¬ 
posed in the NE conditioned process. Our first main 
result is to connect, in Sec. IIIB, the EQ reference pro¬ 
cess and the NE process, and as a consequence, also to 
connect their associated driven processes (see the vertical 
arrows in Table I) . Our second main result is to prove, in 
Sec. IIIC, the asymptotic equivalence between the path 
probabilities of the driven reference process with the NE 
process. This equivalence is at the core of the afore¬ 
mentioned equivalence between the NE micro-canonical 
ensemble and the meta-canonical ensemble. In Sec. IV, 
we comment the structure of the theory starting with a 
short summary in Sec. IV A. We discuss the symmetries 
of the NE potentials and the connection with close-to- 
equilibrium and far-from-equilibrium perturbation the¬ 
ory in Secs. IV B and IV C respectively. We end by illus¬ 
trating our work on a two-level model in Sec.V. 


For the sake of simplicity, we focus on systems exchang¬ 
ing only energy with heat reservoirs. The generalization 
of our results to include matter, volume or other exten¬ 
sive variable exchanges with reservoirs is straightforward 
[36]. 



NE system T2 7^ Ti 

Figure 1. System with M — Q states connected to x = 2 
heat reservoirs at the same temperature Ti = T 2 for the EQ 
reference process, or at different temperature Ti 7 ^ T 2 for the 
NE process. 



Figure 2. Energy lanscape for the x y transition. The 
discrete states x and y represent the locations of the minima 
in the energy landscape. Changing the dilatation factor l\ 
modify the height of all energy barriers for the EQ reference 
process. 


II. EQUILIBRIUM FLUCTUATIONS 
A. Definition of the EQ reference process 

We consider an EQ reference process corresponding to 
a physical system modeled by a continuous-time Markov 
chain with a finite number M of discrete states. This 
system exchanges energy with x heat reservoirs labeled 
by = 1 ... X at inverse temperatures /3i = l/^ksTi), 
with /cs = 1 the Boltzmann constant, see Fig. 1. The 
reference process is at EQ, i.e. all the heat reservoirs 
share the same inverse temperature /3i. We use several 
heat reservoirs to allow different mechanisms of energy 
exchange. As a result, some rare events with net energy 
flow from one heat reservoir to another may occur. The 
system states are generically denoted x, y and z. The 
state at time r is z(t). A system state trajectory during 
time interval [0, t] is denoted [z]. This trajectory includes 
the state z(t) at all time r G [ 0 ,t] and the label z^(t) of 
the reservoir providing the energy at each change of state 
in the trajectory. 

The energy of state x is e^- The probability per unit 
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time of switching from state y to state x exchanging the 
energy tx — Cy with reservoir v is given by the Arrhenius 
transition rates 


^xy — ixy^ 


( 1 ) 


We have introduced the symmetric matrices 7 *^, whose 
(x, y) element yields the coupling with reservoir v for a 
transition from y to x. The (x, y) element of the symmet¬ 
ric matrix d represents the height of the energy barrier 
that must be crossed when the system switches between 
states y and x, see Fig. 2. The dimensionless parameter 
Zi is a dilatation factor that enables to modify the height 
of the energy barriers {h = 1 implies no dilatation). The 
transition rates defined in Eq. (1) verify for all 1 / the local 
detailed balance relation 


(2) 

^yx 

which ensures that the system will reach EQ [37]. The 
reference probability per unit time of escaping from state 
y, given that energy is exchanged with reservoir v, is 
denoted 


per unit time and is regarded as a time-averaged variable. 
At long time, it converges to the step frequency of the 
walker. The probability of v[z] = v, i.e. that the time- 
averaged number of steps v[z\ takes the value v at long 
time t, is Pt{v) ~ The function I is called a large 

deviation function (LDE). It is non-negative and vanishes 
at V = {v{z\)j^, denoting that the ensemble average value 
is the most likely time-averaged v. Small (respectively 
large) deviations correspond to the time-averaged num¬ 
ber of steps that are close-to (respectively far-from) the 
ensemble average value. These events become exponen¬ 
tially unlikely with increasing time for ergodic systems. 
The convexity of the LDE ensures that a large deviation 
is less likely than a small fluctuation. 

Following, we introduce the empirical time-averaged 
variables used to derive our central results. We name em¬ 
pirical variables those that are defined from experimental 
observations of the system and that usually depend on 
the observed trajectory [z]. First, we define the empirical 
occupation in x by 


Px[z] 


1 

t 





( 4 ) 


= (3) 

x=^y 

such that each column of the matrix sums to zero 
as required for continuous time Markov chains. The ref¬ 
erence transition rate matrix k = k^'''> returns the 

transition probabilities per unit time disregarding the 
reservoir involved in the energy exchanges. Similarly, 
is the total escape-rate vector. As a conven¬ 
tion, we drop the subscripts of vector or matrix elements 
to refer to the whole vector or matrix and use bold face 
letters for matrices. We denote the ensemble average over 
all trajectories [z] generated with dynamics correspond¬ 
ing to k with the brackets (. • ■)fc. 

B. Large deviations of empirical time-averaged 
variables 

Throughout the paper, we assume that the long-time 
statistics of time-extensive variables obey a large devia¬ 
tion principle. For instance, let z(t) be the position at 
time t of a random walker on a one-dimensional circu¬ 
lar lattice and X[z] the number of steps the walker takes 
during the trajectory [z]. We remark that the variable 
A is a functional of the trajectory [z] that is a realiza¬ 
tion of a stochastic process. X is not a random variable 
in itself. When X is not evaluated on a trajectory, it 
refers either to the physical variable “number of step” or 
to a numerical value of this variable. The variable X[z] 
is time-extensive since X[z] -\- X[z'] = X[z,z'], where 
[z, z'] denote the trajectory made with [z] followed by [z']. 
Then, the number of transitions typically increases with 
time. Accordingly, v[z\ = X[z\/t is the number of steps 


where 5 is the Kronecker symbol. Given the probability 
of each state being gathered into the column vector p = 
(pi, • • • jPmY, the Shannon entropy s = s{p) is 

sip) =-'^Pxlnpx =-{p'^-Inp), (5) 

X 

and the energy e = e{p) is 

e(p) = ^ exPx = -P, (6) 

X 

with the central dot denoting the matrix product and f 
the transposition. The time-averaged energy along tra¬ 
jectory [z] can be written e[z] = e(p[z]), and similarly for 
the entropy. Second, we define the empirical transition 
probability from y to x induced by reservoir v 

^ ' ^x,z(T+dT)^y,z(T)^i',i'{T)^ ij) 

rG[0,t] 

where the sum is over all time r at which the system 
changes from state z(t) to state z(r -I- dr), exchanging 
energy with reservoir v{t). Given a transition probability 
ujxy for each possible transitions, the current of energy 
received from reservoir v by the system is 

j'^i‘^)=\'^{^xy-‘^yx)i<^x-ey)- ( 8 ) 

x,y 

Its empirical value during trajectory [z] is written jp[z] = 
jp{uj[z]). These time-averaged currents describe the anti¬ 
symmetric part of fluctuations since they change sign 
upon time-reversal of the trajectories. On the contrary, 
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the weighted frequency of interaction with reservoir z/, 
named activity for short and written 

f^{u}) = - ^ + ujy^) dxy, (9) 

x,y^x 


describes the symmetric part of fluctuations. Indeed, 
fv[z] = f^{u}[z]) does not change sign upon time-reversal 
of the trajectory [z]. When the activity is low (high), the 
system either changes of state less (more) frequently or 
mostly switches between states with low (high) d^y. The 
term “activity” was proposed to qualify the symmetric 
part of the fluctuations in [22, 38-41], see also references 
therein. Let us finally remark that, in the definitions of 
the energy currents and activities, the one-half factor is 
just a symmetry factor since we can sum over transitions 
disregarding their directions y) only one di¬ 

rection (X)a;>y)- Half of the first sum is equivalent to the 
second sum. 


C. LDF of energy currents, activities and 
occupation from a variational approach 


At long time t, the probability of observing an em¬ 
pirical transition probability <jj{z\ = (jJ and an empirical 
occupation p[z\ = p is 

Pt{u>,p) ~ ( 10 ) 


From the work of Maes and Netocny [27], Wynants [42] or 
Bertini et al. [43], the LDF I{uj,p) of the empirical tran¬ 
sition probabilities and occupations for the continuous¬ 
time Markov chain with generator k is 


I{(^,p) = 


x,y^x,i> ^ 




Py 


xy 


■ UJ In -—— , 

^y 71 

^xyPy J 


( 11 ) 


where the sum is over v from 1 to y and all couples 
{x, y) such that x ^ y. The derivation of Eqs. (10-11) is 
reproduced in Appendix A. 

In Ref. [27], the LDF of the occupation and probabil¬ 
ity current was obtained from a constrained optimiza¬ 
tion problem constructed with I{u},p). This procedure, 
called “contraction” [4] , is equivalent at the level of prob¬ 
abilities to marginalize Pt{<jJ,p) to obtain the probability 
of currents and occupations. We now proceed to the 
contraction of I{u},p) to obtain the LDF of energy cur¬ 
rents, activities and occupations denoted L{j,f,p). The 
long-time asymptotic approximation of the probability 
Pt{j, f,p) that j[z] = j, f[z] = f and p[z] = p at time t 
will then read 

PtijJ,p)^c^ ( 12 ) 

t—¥oo 

We prove in Appendix B a sharper approximation of this 
probability that involves a pre-exponential factor dictat¬ 
ing the thermodynamic behavior: it leads to the statistics 


of the usual EQ thermodynamic variables that only de¬ 
pend on the system state, such as energy for instance. 
It is the first correction to the exponential decay of non¬ 
typical time-extensive variables after a long time. This 
prefactor was first obtained in Ref. [41], but we pro¬ 
vide in Appendix B a logically independent derivation 
(though restricted to the large time limit) that involves 
some results of Sec. III. 

At long time, the energy current j, activity / and oc¬ 
cupation p mainly appear thanks to the most likely event 
producing them. The probability of this event is associ¬ 
ated with smaller values of I{uj,p) with u} constrained 
by the value of the energy currents and activity. For 

this reason, and in virtue of the contraction principle, we 
minimize I{lj,p) under the energy currents constraint 

>=>(w), (13) 

for 1 / > 1, because current conservation imposes ji = 
We also impose the activity constraint 

fu = (14) 

In addition, the probability currents should be compati¬ 
ble with the conservation of the norm of the occupation 
vector, i.e. for all y 

i:k,-o=o. (15) 

x,v 

To perform our optimization problem, we use the follow¬ 
ing cost function 


P{u},p) = I{u:,p) +^a^[j^ - j^{u})] 

V 

I/' x,y,i' 

where a,y, b^, and Uy are Lagrange multipliers that will 
be chosen to satisfy the constraints of Eqs. (13-15). We 
choose oi = 0 so as not to constrain the current ji that 
is already set by the current conservation law. We now 
minimize the function P with respect to lj, calculating 
dPjduj^y = 0 to get 

0 = In — ai/icx Cy) biidxy -t- (uy ^x), (1"^) 

^xyPy 

where we have used Eq. (11) and Eqs. (13-14). There¬ 
fore, the optimal transition probability in terms of the 
Lagrange multipliers satisfies 

^xy = KyPy, ( 18 ) 

where we have introduced K" = K''{a,b,u), the transi¬ 
tion probability for mechanism v divided by the empirical 
occupation of the state before transition. Its off-diagonal 
elements are 


ATI 


^1/ ^a^{€^-ey)+b^d^y+U:^ 
xy 


(19) 
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or more explicitly using Eq. 1 

^-(l3il2-a^)(e^-ey)-{fiih-b^)d^y+u^-Uy 

and the diagonal elements are 

= (21) 
x^y 

such that any column of any matrix K'' sums to zero. We 
remark that the matrices K’' satisfy a modified detailed- 
balance relation 

K'' 

In —^ = {2ay, - /3i)(ex - e^) -f 2{ux - Uy). (22) 

-^yx 

In this local detailed balance, the Lagrange multiplier ay, 
biases the inverse temperatures /3i to make typical the 
energy exchanges corresponding to the energy currents 
constraint. The reservoir v behaves as if it had the tem¬ 
perature = /3i — 2ai, in order to satisfy the current 
constraint. Thus, the variable 

2a, = /3i - /3, (23) 

is an affinity [44-46], also called thermodynamic force 
[47, 48]. Notice that oi = 0 as required. The similarity 
between Eqs. (1) and (20) indicates that we can also in¬ 
troduce new dilatation factors Z, such that the dynamical 
bias 

bu = Pih — Pv^v (24) 

gives the modification of the dynamics in order to satisfy 
the activity constraint. Finally, we call the variable u the 
drift because it acts like a force biasing each transition. 

The explicit solution w of our variational problem 
dJ^ = 0 is now almost reached. The next step is to use 
the constraints of Eqs. (13-15) to obtain the Lagrange 
multipliers. More explicitly the constraint equations are 

jv — 2 E! i^xyPy ~ ^yxPx) i^x ~ ^y), (25) 

x,y 

fv = 2 E! i^xyPy + ^yxPx) dxy, (26) 

^,V 

0 = K-p, (27) 

where K = K'' is the generator of the driven refer¬ 

ence process [29, 31]. For the third equation, the con¬ 
servation law of the probability current of Eq. (15) is 
reformulated as a requirement that the empirical occupa¬ 
tion p is the stationary probability of the continuous-time 
Markov chain with rate matrix K = K{a,b,u). Invert¬ 
ing these three equations gives the vectors a, b and u as 
a function of (j, /, p). 

The final step to obtain the asymptotic probability of 
energy currents, activities and occupations is to write the 
LDF of Eq. (11) at the optimal transition probability of 
Eq. (18). This leads to 

L{j, f,P) = ■ j ■ f {X- A)^ ■ p, (28) 


where we have used the anti-symmetry of — Cy or sym¬ 
metry of dxy in the exchange of x and y to make explicit 
the dependence in j and /. We also used Eq. (15) to get 
rid of the term involving Ux — Uy 

D. Scaled CGF of energy currents, activities and 
occupations from a variational approach 

In the previous section, we have obtained the LDF L 
from the solution of an optimization problem. From now 
on, and for the remainder of the paper, we assume the 
convexity of the LDF. Our aim here is to derive the scaled 
CGF conjugated to L from a variational approach, using 
the fact that LDF and scaled CGF are conjugated by 
Legendre transformation [4, 5, 31]. On the way, we obtain 
useful properties associated to the canonical structure. 

The scaled CGF of the energy current, activity and 
occupation is defined by 

F(a', 6 ',m')= lim 1 In 

(29) 

and is the Legendre transformation of L 

F(a', b', m') = max [a'''' ■ j + b'^ ■ f + ■ p - L{p,j, /)] . 

pJJ 

(30) 

The maximum on j and / is reached for a' = a and b' = b, 
and the scaled CGF becomes 

r(a, 6 , m')= max [(m'-I-A — A)^ • p] , (31) 

p I K p=o ■' 

where the maximum is taken over all occupations with 
the Lagrange multiplier u in the generators of the driven 
process K tuned such that K ■ p = 0. An alternative 
expression of the scaled CGF of energy current, activity 
and occupation is 

r(a, &, to') = max [(to'-I-A — A)^ • p] (32) 

with p the stationary probability associated to K. From 
the optimal drift u = u{a, b, to') realizing the maxi¬ 
mum in Eq. (32), we introduce the escape weight to = 
m(a,b,u) giving the value of to' for given {a,b,u). In 
Eq.(G14) of the appendix of reference [28], Nemoto and 
Sasa gave the scaled CGF of energy current from a vari¬ 
ational expression analogous to our Eq. (32). We recover 
their result taking 6 = 0 and to' = 0. We further com¬ 
ment Eq. (32) noticing that the maximum is reached for 
u satisfying 

F(a, 6 , to) = ITly Ay — Ay, (33) 

for all y. This equation allows us to derive the following 
escape-rate rule 

TOy Ay Ay - TO-j; A^, A 3 ;, (^^) 

that can be related to the exit rate constraint of Refs. [24, 
25, 29] taking to = 0. To prove Eq. (33), we introduce 


7 


the tilted operator k = K{a,b,m) for the EQ reference 
process 

i^vv = -^Ky + my, (35) 

x^y,v 

K^y = . (36) 

V 

The generator of the driven reference process iC is con¬ 
nected to this tilted operator by 

d^xy — ^ ^ rixyG y (^y “f ^y '^y) ^xy- (^^) 

Using this equation and K ■ p = 0, we find 

^ ^ e ^ f^xyG ^Py — (jGlx “t“ d^x ^x') Pxi (^^) 

y 

after summing over x and maximizing over rt, it follows 
from Eq. (32) that the drift giving the maximum satisfies 

=r. (39) 

x,y 

By definition [49], the scaled CGF r(a, b, m) is the high¬ 
est eigenvalue of n. Then tt^, = ±e““’/Z(M) is the nor¬ 
malized left eigenvector of k with Z{u) a normalization 
constant such that J2x'^x = 1- The vector r = -p is 
a right eigenvector with iTxy = T^xdxy Its norm is set by 
'^x '^x'l'x = J2xPx = 1- Notice that we cannot determine 
from the values of u the sign of each component of the 
vectors tt and r, but their x components share the same 
sign. Now, summing Eq. (37) over x leads to Eq. (33) 
since Yjx ^xy = 0 and = T. 

Then, the optimal drift u = u{a,b,m), leading to 
the maximum in Eq. (32), is simply obtained from the 
left eigenvector of the tilted operator by ln|7ra;| = Ux — 
ln|Z(u)| up to a constant that plays no role, since only 
differences of drifts matter. The drift makes the escape- 
rate rule holds true and, using Eqs. (28) and (33), leads 
to the Legendre structure that one expects for LDFs and 
scaled CGFs. Finally, from Eqs. (33) and (37), we recover 
the results of Refs. [29-31, 34] in which the generator K 
of the driven process corresponds to the Doob’s transfor¬ 
mation of the tilted operator n 

^xy — I'^a: l^a:y ["^y 1 ^dxy (^H) 

Notice that in Refs. [29-31], the right eigenvector of the 
tilted operator is used in the Doob’s transformation in¬ 
stead of the left one, since the tilted operator in these 
references is the adjoint of n. 


III. NON-EQUILIBRIUM FLUCTUATIONS 


reservoir v at different temperatures = Pi — 2ai, and 
with different dilatation factors G related to dynamical 
bias by = PPi — The elements of the rate ma¬ 
trices are 

(41) 

Accordingly, the escape rate Xy = Xyia^^b^) from state 

y is 


K^T.~kxy = -Ky ( 42 ) 

x=Xv 

We define a total rate matrix by fc = Yv ^ total 

escape-rate vector by A = Yu ■ These rates are func¬ 
tions of the affinities and dynamical bias; their analogs 
for the reference process are recovered at the point of 
vanishing of a and b, namely k = fc(0,0) and A = A(0,0). 
For the NE process, the state at time r is z{t). A system 
state trajectory during time interval [ 0 ,t] is denoted [z]. 
The ensemble average over all trajectories [z] generated 
with dynamics corresponding to fc is (.. .)^. 


B. Mapping typical NE fluctuations on rare EQ 
fluctuations 


We now connect the energy currents, activities and oc¬ 
cupations statistics for the EQ process with the statistics 
of the same variables for the stationary NE process. This 
mapping involves the escape-rate change c = c(a, b) de¬ 
fined by 


c = A - A, (43) 

that is zero at vanishing affinities and dynamical biases. 
We emphasize that c cannot be adjusted independently of 
a and b. This means that the affinity and the dynamical 
bias are the central variables in determining the NESS 
reached by the system. 

To connect EQ and NE fluctuations, one needs to 
redo all the calculations of sections 11C and IID, but 
for the NE process, introducing the NE scaled CGF f = 
f (a, b, fh) and LDF L — L{j, f,p), the NE tilted operator 
k = k{d,b,m), the generator of the NE driven process 
K = K{d, 6 , u) and associated escape rate A = A(a, b, u), 
the affinity increment 2 a, the dynamical bias increment 6 , 
the NE drift u and the NE escape weight fh = fri{d, b, u), 
all denoted with a bar to distinguish them from their 
equivalent for the EQ reference process. One obtains all 
these objects replacing fc by fc and the Lagrange multipli¬ 
ers (a, b, u) by (a, b, u) in all the definitions. For instance, 
for the NE tilted operator, we have 


A. Deflnition of the NE process 

The NE process is defined by the rate matrices k’' = 
k''{av,bu) associated to energy exchanges with each 


^vv ~ ^ ( ^xy T 'eriy, 

(44) 

x^y,u 


uv „a^{e-^-^y)+bud:^y 
i^xy — / . f^xyG 

1/ 

(45) 


Notice that we call 2a an affinity “increment” since we 
already deal with a NE process: a deviation from the 
typical current is associated with an “increase” of affin¬ 
ity that will make this fluctuation typical. For the same 
reason, the dynamical bias b is also qualihed as an incre¬ 
ment. 

The mapping between EQ and NE fluctuations now 
comes from the connection between the EQ and NE tilted 
operators 

K{d, b, fh) = K{d -I- a, 6 -I- 6, m -I- c), (46) 

that we obtain by comparing Eqs. (35-36) with Eqs. (44- 
45). Hence, the same symmetry exists between the eigen¬ 
values and between the eigenvectors: the full spectrum of 
the two operators is connected. In particular, the scaled 
CGFs are connected by 

f (a, 6 , m) = r(d -I- a, 5 -I- 6 , m -I- c), (47) 

and, from the Legendre transformation, the LDFs verify 
f,p) = L{j, f,p)-a^-j-b^-f-c^-p. (48) 

The left eigenvectors of the tilted operators satisfy 

7r(d, 6, fh) = 7r(a + a,b + b,fh + c) (49) 

or equivalently 

u(a, b, fh) = u(a + a,b + b^fh + c). (50) 

The mapping also holds for the right eigenvectors and 
this leads to 

p{a,b,fh)=p{a + a,b + b,fh + c). (51) 

Finally, the generators of the driven processes also verify 

K{a,b,u) = K{a + a,b + b,u), (52) 

where u and u are respectively the left- and right-hand 
sides of Eq. (50). 

Thus, the EQ and NE processes are tightly connected 
and one can focus on the EQ process’ fluctuations only: 
Eq. (47) shows that the statistics of energy currents, ac¬ 
tivities and occupations for any NE process with affinity 
2a and dynamical bias b is known from the statistics of 
the same variables computed for the EQ process. Indeed, 
the derivatives of Eq. (47) with respect to a,b or fh eval¬ 
uated in (a, b, fh) = (0, 0, 0) yields the NE cumulants of 
the energy currents, activities and occupations from the 
scaled CGF for the EQ reference process, e.g. for we 
have 

(>[^])s = ^(0,0,0) = ^(a,5,c). (53) 

Notice that evaluating Eq. (47) at the point of vanishing 
of (a, b, fh) returns by definition of a scaled CGF 

0 = f(0,0,0) = r(a,6,c), (54) 


for all a and b, with c = (A — A). Accordingly, the total 
derivatives of r(a, &, c) with respect to a or 6 also vanish 
exactly such that T remains constant and equal to zero in 
the direction (a, 5, c). We call the subspace where T van¬ 
ishes the physical system subspace: each point (a, &, c) in 
this subspace defines a precise physical process with affin¬ 
ity 2a and dynamical bias b. The function T includes the 
full thermodynamic information on any system defined 
with the same energy levels e, coupling matrices 7 “^ and 
energy barriers d (up to a reservoir specific dilatation), 
and so does the LDF L. One simply changes the degree of 
NE or the type of dynamics, encoded into the dilatation 
factors, by moving into the physical system subspace. 

We end by remarking that the idea of mapping EQ 
and NE fluctuations was first proposed by Andrieux in 
Refs. [50], but for the statistics of energy currents only. 
However, this mapping had no concrete application since 
the NE statistics of the currents were needed to define 
the EQ dynamics involved in the mapping. On the con¬ 
trary, the mapping of Eq. (47) and (48) is explicit, with 
the price that, when comparing with Ref. [50], the EQ 
statistics of activities and occupations must be known in 
addition to the energy currents statistics. 

C. Asymptotic equivalence of the driven reference 
process and the NE process 

We now discuss the asymptotic equivalence of the 
driven reference process and the NE process. We first 
prove the equality of their escape rates and on the way 
give a slightly simplified expression of L. Using this re¬ 
sult, we demonstrate the equivalence of the path proba¬ 
bilities of the driven reference process and the NE pro¬ 
cess. 

From Eqs. (33) and (54), we find c -I- A — A = 0. This 
leads with Eq. (43) to the equality of the escape rates of 
the driven reference process and the NE process 

A = A, (55) 

even though these two processes are different in general 
due to the drift u, i.e. K^y k^y ii x ^ y. As a conse¬ 
quence, the LDF is written as 

L{jJ,p) = a^-j + b^ ■ f + c^-P- (56) 

The equality of the escape rates indicates that the driven 
reference process and NE process look alike. Their gen¬ 
erators are connected by the similarity transformation 

iT = [ttI • fc • |7r|"\ (57) 

that follows from the comparison of Eqs. (20) and (41). 
We denote | 7 r| the positive and diagonal matrix obtained 
by taking the absolute value of the elements of tt. The 
equality of the diagonal part of the Markov matrices 
of the two processes is granted by Eq.(55). From this 
similarity transformation, one can show the asymptotic 
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equality of the path probabilities associated to each pro¬ 
cess 


VkIv]^^ r-k[y], (58) 

for any trajectory [y]. We have defined the path proba¬ 
bilities knowing the initial state y{0) 

rk[y]=exp(- f dTXy(^r)) H Kirldr)y(ry (59) 

^ do / 

for the NE process and 

VK[y] = exp(- f drAy^r)] H K{r+dr)y(Ty ( 60 ) 

^ do / 


for the driven reference process. In these equations, 
the product applies for all times r at which the system 
changes of state during the trajectory [y], with y(r) (re¬ 
spectively y{T+dT)) the system state before (respectively 
after) the transition at time r. The exponential terms ap¬ 
pearing in these two path probabilities are equal. Con¬ 
cerning the product terms, they differ from boundary 
terms only 


J_ ^y{T+dT)y(T) J. l^y(r-t-(iT)y(r) 


1-1 


r I j 


~ Ky(t)l l^n \(T+dr)y(T) J l^y(0)l 

(61) 


Then, the path probabilities of the driven reference pro¬ 
cess and NE process verify 


lim - In 

t—^OO t 


niy] 

VK[y] 


= 0 , 


(62) 


and are asymptotically equivalent [29]. Since the driven 
reference process is the dynamics that typically repro¬ 
duces the conditioned reference process, we conclude that 
there is an ensemble equivalence between the NE process 
and the conditioned reference process. This central re¬ 
sult is similar to the path-ensemble equivalence derived 
in Refs. [29, 30]. In Appendix C, we show that the NE 
process is asymptotically equivalent to the canonical pro¬ 
cess that is defined by exponentially weighting each tra¬ 
jectory, even though these two processes are not exactly 
identical. 


IV. DISCUSSION AND GENERAL SUMMARY 

In Sec. II, we have studied the fluctuations of an EQ 
system exchanging energy with several heat reservoirs at 
the same temperature. We have seen that energy may 
spontaneously flow from one reservoir to another, even 
if it does not on average. Each of these current fluctu¬ 
ations has been associated to a temperature difference 


that would typically reproduce it. Similarly, we have 
shown that a fluctuation of the activity of the exchanges 
with each reservoir would be typically reproduced by di¬ 
lating the appropriated energy barriers. From these ob¬ 
servations, we have identified two couples of conjugated 
variables and provided the corresponding LDF and CGF 
from a variational approach. 

In Sec. Ill, we have considered the fluctuations of the 
system defined in Sec. II, but driven out-of-equilibrium 
by temperature differences between the heat reservoirs. 
We have found an exact mapping between the statistics 
of the energy currents, activities and occupations for the 
EQ and NE systems. We have also discussed the asymp¬ 
totic equivalence of the trajectory ensembles generated 
by the conditioned EQ process and the NE process. From 
the existence of the mapping between EQ and NE sys¬ 
tems, we have concluded that the study of a NE system 
amounts to the calculation of the probability of rare fluc¬ 
tuations of the same system at EQ. Now that the distinc¬ 
tion between the dynamical fluctuations of EQ and NE 
systems has been dispelled, we come back to the results 
of Sec. II and summarize the canonical structure satisfied 
by the two ensembles of NE systems. 


A. Summary of the NE canonical structnre 

The ensemble of systems in contact with several heat 
reservoirs at different temperatures is called the meta- 
canonical ensemble. The trajectories of the systems in 
the meta-canonical ensemble are generated by the NE 
process with generator k. All the systems in this ensem¬ 
ble have the same energy levels ex, and the same dynam¬ 
ical parameters, i.e. energy barriers dxy and couplings 
with the heat reservoirs "fxy convention, the heat 
reservoir of smallest temperature is the reference reser¬ 
voir (jz = 1) such that all the affinities 2ajy = Pi — Pu 
are positive. Notice that the temperature of the refer¬ 
ence reservoir sets the energy scale and has no physical 
relevance. On the opposite, the affinities ai, are the cen¬ 
tral variables of the meta-canonical ensemble that are set 
by the environmental constraints. The affinities are nat¬ 
urally conjugated to the energy currents. However, we 
know from the previous sections that considering (o, j) 
as the unique couple of conjugated variables does not 
afford to study all NE systems from the same NE po¬ 
tential. Intuitively, a change of an affinity also impacts 
the system activity and the occupation of the various 
states. Hence, we have introduced additional intensive 
variables to take into account these effects separately: 
the dynamical biases connected to the dilatation factors 
of the energy barriers and the escape weights modifying 
the escape probability of each state. These two intensive 
variables cannot be adjusted independently of the affini¬ 
ties if we want to avoid a change of the system dynam¬ 
ics: no dilatation should be applied to the energy barriers 
(?y = 1 for all v) yielding to dynamical biases that are 
equal to the affinities (6 = 2a); the dynamics should con- 
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Table II. EQ and NE thermodynamic potentials. 


Ensemble 

Micro-canonical 

Canonical 

NE micro-canonical 

Meta-canonical 

Potential 

s =-YjxPx\ripx p = 

-In {exp{-Pie[z]))^ 


r(a, b, m) 

Variational principle 

max 

min 

max 

max 

Free variables 

Pi 

e 

a, b, m 

j, f, p 

Constrained variables 

e 

pi 

j, f ,P 

a, b, m 

Physical system subspace 

— 

— 

m{j, f,p) = c(a(j, f,p), h{j, f,p)) 

m = c(a, b) 

No dilatation space 

— 

— 

b{jj,p) = 2a{j,f,p) 

b — 2a 

Legendre structure 

s -b v? = 

Pie 

L-\-T = a'-j-\-b'-f-\- m' • p 


serve the norm of the occupation vector imposing that an 
affinity must be associated with an escape weight equal 
to the escape-rate change m = c(a,2a). Therefore, in 
the meta-canonical ensemble, the environment sets the 
affinity vector a which in turn constrains the dynamical 
intensive variables, namely the dynamical biases and the 
escape weights. The NE potential of the meta-canonical 
ensemble is the CGF of energy currents, activities and 
occupations r(a, 6, m). It vanishes for all a when b = 2a 
and m = c(a, 2a), but its partial derivatives with respect 
to a, 5 and m produces all the NESS cumulants of energy 
currents, activities and occupations for any affinity. For 
instance, the thermodynamic behavior follows from the 
NE equations of state 


ar 

dai, 



dr 

db,y 

(9r 

drrix 


= fu, 

,m 

= Px, 


(63) 

(64) 

(65) 


where the subscripts on the vertical bars indicate vari¬ 
ables that remain constant when taking the partial 
derivative. We denote the vector a without the I'th 
component. The cumulants of EQ thermodynamic vari¬ 
ables are obtained with the NESS occupations defined 
by P* = p{a, 2a, c) that only depend on x — 1 affinities. 
The mean energy in the NESS is {e[z])j;, = e(p*), and the 
mean entropy is (s[z])fc = s{p*). 

The ensemble of systems conditioned on the energy 
currents they received from their environment is called 
the NE micro-canonical ensemble. The trajectories of 
the systems in this ensemble are generated by the EQ 
reference process with generator k filtrated to achieve 
the condition on the energy currents. The physical im¬ 
plementation of systems in the NE micro-canonical en¬ 
semble would require the existence of energy sources with 
no fluctuations. These sources will very likely not exist 
in practice [51], even though this problem is not specific 
to NE ensembles (see for instance page 83 of Ref. [3] for 
an example in EQ thermodynamic theory). If we assume 
that an energy current can be imposed from the out¬ 
side, the activities and the occupations must take pre¬ 
cise values so that the system can sustain the energy 
current. On the opposite, the conjugated intensive vari¬ 


ables become free to fluctuate. The relationship between 
currents, activities and occupations is obtained from the 
correspondence between the conjugated variables {j, f,p) 
and {a,b,m), as summarized in Table II. The NE micro- 
canonical potential is the LDF L{j, f, p) and the statistics 
of the intensive variables (a, b, m) follows from its partial 
derivative 


dL 

dju 

dL 

dfu 

dL 

dpx 



— Cti/, 

(66) 





— 

(67) 





= rUx- 

(68) 





We proved in sections II and IIIC the equivalence of 
the ensembles of trajectories generated by the NE process 
and the conditioned EQ reference process assuming that 
the NE potentials are convex. Accordingly, the meta- 
canonical ensemble and NE micro-canonical ensembles 
are ensemble equivalent. In other words, systems sub¬ 
mitted to temperature gradients are equivalent, at the 
thermodynamic level, to systems subjected to stationary 
energy injection (and extraction). By construction, the 
NE potentials are conjugated by Legendre transforma¬ 
tion 


L{j, f,p) + r(a, b,m) = ■ j -\-b^ ■ f -\-m^ -p, (69) 

and the NE stationary state can be obtained from a vari¬ 
ational approach. If we consider • j -|- 6^ • / -|- ■ p — 

r(a, b, m) as the potential L that would be obtained from 
Eq. (69) by assuming the independence of the conjugated 
variables, then the NESS affinity, dynamical bias and es¬ 
cape weight reached by the system at constant imposed 
energy current j, activity /, and occupation p maximize 
this potential in the subspace of constant {j, f,p): 

(a, b, m) = argmax [a^ • j -I- 6^ • / -I- ■ p — T{a,b, m)] 

(70) 

which are exactly Eqs. (63-65). The same argument holds 
the other way around. If we consider ■j-\-b^ ■ -p— 

L{j, f,p) as the potential T that would be obtained from 
Eq. (69) assuming the independence of the conjugated 
variables, then the NESS energy currents, activities and 
occupations reached by the system at constant imposed 
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affinity a, dynamical bias b, and escape weight m maxi¬ 
mize this potential in the subspace of constant {a,b,m): 

(j, f, P) = argmax [a^ ■ j + b'^ ■ f + m'^ ■ p - L{j, f, p)] 

(71) 

which are exactly Eqs. (66-68). 


B. Symmetries of the NE potentials 

The metacanonical potential is even under the sign 
change of all affinities. We prove in Appendix D that this 
symmetry leads to the fluctuation theorem (FT), a fun¬ 
damental result regarding the asymptotic statistics of en¬ 
tropy production first studied in Refs. [52-54]. Another 
fundamental symmetry is obtained from the equality of 
second derivatives of the NE potentials. This symmetry 
is the NE equivalent of the Maxwell relations and reads 
as 


d^T 

dhadh'^, 


d^L 

dh'^,dK dv^dv'^, 


d^L 

dv'^,dva 


(72) 


where h and h' are two vectors in (a, 5, m) and similarly 
V and v' in (j, f,p). The subscripts a and a' indicate two 
arbitrary components of these vectors. At EQ, Maxwell’s 
relations deeply constrain the number of EQ response 
coefficients that should be introduced to completely de¬ 
scribe a system. Here, they constrain the derivatives 
of the non-linear functions giving, for instance, the cur¬ 
rents in terms of the affinities. In the close-to-EQ limit, 
Eq. (72) implies that the linear response matrix is sym¬ 
metric, or in other words it implies the Onsager reci¬ 
procity relations [55, 56], as we will see in the next sec¬ 
tion. 


C. NE linear response theory 

We study the linear response of a system in an arbi¬ 
trary NESS and further perturbed by a change of tem¬ 
perature /3i, —>■ /3(, = /3i, -I- A/3y or of dilatation factor 
Iv ^ = lu + Al^. More precisely, we want to de¬ 

termine the change of the energy currents and activities 
when the half affinities Qi, = (/3i — /3^)/2 and dynamical 
biases b,^ = are slightly changed to the new 

values -I- Aa^ and b^, + Ab^. We assume that li and 
/3i do not change during the perturbation. Then, the 
perturbations are written as 

Aa, = -A/3,/2, (73) 

A6, = (—/3, -l- 2Aa,) AQ -I- 2Z,Aa, 

— —/3i/AZ, -|- 21,Aa,, (1'4) 

at linear order. We remark that the dynamical biases 
change when perturbing the affinities, but the converse 
is not true. 


A Taylor expansion of the meta-canonical potential T 
gives the following quadratic function 


r(a -I- Aa, b + Ab, m + Am) = 

r{a,b,m) + Ah^-VfiT 

h—a,b,m 

+ ^ ^ Ah^-Vhh'T-Ah', (75) 

h —a,b,m 
h'—a,b,m 


where Am is not yet specihed. We have used the short 
notations for the derivatives of the meta-canonical po¬ 
tential 

dT 

(V/,r)^ = — (a,6,m), (76) 

The perturbation induces a variation Aj of the energy 
currents. A/ of the activities or Ap of the occupation. 
Taking the partial derivative of Eq. (75) with respect 
to Aa, A6, or Am and evaluated in Am = Ac, with 
Ac the variation of the escape-rate change due to the 
perturbation, leads to the linear response equation 


(Aj\ 

' VaaE 

VafeE 

VamE ■ 

A/ U 

VfioE 

Vf,bE 

VbmE 

\ApJ 

_ V^aE 

Vmf)E 

V r 

v mm^ 


Aa \ 

Ab . (78) 
Ac / 


From Eq. (72), the response matrix above is symmetric 
even close to an arbitrary NESS. However, the chain rule 
yields 


Ac = VaC • Aa -|- VfeC • Ab, (79) 

and the variation of the currents and activities becomes 


Aj = (Vaah -t- VamF • VaC) • Aa 

+ (VafcT -|- VomF • Vfcc) • Ab. (80) 
Af = (VabT -I- VbmF • Voc) • Aa 

+ (VbbF -|- VbmE • Vf,c) • Ab. (81) 

The response matrix defined from Eqs. (80-81) is no 
longer symmetric in general as already emphasized in for¬ 
mer works on NE linear response theory [38, 39, 57-62]. 
The second derivatives of the meta-canonical potential 
appearing in Eq. (80) are 

(Vaah),,, = lim t {(/,[%,'[z])fe - (j,[z])fe (j,'[2])fc}, 

t—yoo 

(VahE)^^, = lim i{(>[z]/,.'[2])fc - {j4z]}k {U'[z])k} . 

t—yoo 

(Va„^^)^,^ = lim t {(j,[zK[z])fc - {ju[z])k {Pa:[z])k} > 

t—yoo 

(82) 

and correspond respectively to the current-current, the 
current-activity and the current-occupation covariances 
in the unperturbed NESS [63]. In addition to the above 
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covariances, the response functions include another con¬ 
tribution involving the derivatives of the escape-rate 
change c. Since the escape-rate change satisfies 

r)r 

y 

the unperturbed mean-occupation multiplied by this 
derivative returns the unperturbed mean energy current 


dc. 


^ Kx ^M)k (ey-ex) = (>[^])fc ■ 




(84) 


Therefore, the response to the affinity perturbation is 


that the energy currents and occupations are decoupled. 
The Yamamoto-Zwanzig formula follows from Eq. (80) 

Aj = ^ • (/?! - /3'), (89) 

where Vooh is given in the first line of Eq. (82) with 
EQ averages instead of the NE averages • •)k- 

Therefore, we recover the Onsager reciprocity relations 
from the NE Maxwell-relations. 


V. ILLUSTRATIVE EXAMPLE: A TWO-LEVEL 
SYSTEM 


(Vaar -|- Va^r ' aC)uu' 

= /-I" ^ • ^)) j ■ 

(85) 

As expected, the response has an additive structure with 
an equilibriumlike part given by a currents correlation 
function, and a NE part corresponding to a current and 
traffic-excess correlation function. We call traffic-excess 
the derivative of the empirical escape rate [z] ■ X with 
respect to the perturbed variable [38, 39]. Similarly, the 
response of the energy current to a perturbation of the 
dynamical bias in the second line of Eq. (80) has two 
parts with an activity-current correlation function and a 
current-traffic excess correlation function. 

As regards the perturbation of an EQ system, i.e. all 
are close to the reference inverse temperature /3i, one 
recovers the Yamamoto-Zwanzig formula expressing the 
response coefficients to a temperature perturbation from 
the covariances of energy currents [64, 65]. In order to 
see this, let us Hrst consider a reference system at EQ 
only perturbed by a change of the dilatation factors, i.e. 
Aa = 0 and Ab = —/3iAZ. Thanks to Eq. (80), the 
variation of the energy currents is written as 

Aj = (VabE -I- VamE • V^c) • A5 = 0. (86) 


We now illustrate our results on a two-level system 
with states z = 1,2 and mechanisms ^ = 1, 2, • • • , y en¬ 
abling energy exchanges with y different heat reservoirs. 
The coupling strength with reservoir v is denoted in 
this section since it is not a matrix but a vector when 
there are only two states. The energy states are ei and 
62 - Let e± = ei ± €2 to shorten notations. The tran¬ 
sition rate matrix of the EQ reference process for each 
mechanism v is 


k‘' 


01 ^- 01 ^- 
-IvG = 7^/6 2 

file- _ fiie- 

. 2 2 


(90) 


where we assume vanishing dilatation factors h, see 
Eq. (1). The rate matrices for the NE system are 

p{Pl/2-a^)e-+b^t+ -{l3i/2-a^)e-+b^e+ 

^ ~ .y^g(/3l/2-a„)e_+&,,e+ _^^p-(l3i/2-a^)e-+b^e+ J 

(91) 

if we chose di 2 = e+. The escape-rate changes for this 
model are 

Cl = ^ (I - , (92) 

V 

C2=Y, (1 - . (93) 


It vanishes for any perturbations A6 since no mean en¬ 
ergy current exists at EQ. Thus, we Hnd 


The tilted operator k = K{a,b,m) for the EQ reference 
system is 


VoftE -f VomE • Vf,c — 0, (87) 


if the derivatives are taken in a = 0. This removes the 
contribution due to the dynamical bias from the EQ re¬ 
sponse. Another contribution disappears in the close-to- 
equilibrium limit due to the decoupling between occupa¬ 
tions and energy currents [42]. Indeed, from the symme¬ 
try of the meta-canonical potential with sign change of 
the affinities, namely E(a, 6, m) = E(—a, 6, to), we have 


a^E 

dttpdmx 


(a, b, to) 


a^E 

dapdrux 


(—a, b, to). 


( 88 ) 


Accordingly, VamE = 0 if the derivatives are taken in 
a = 0. From the third line of Eq. (82), we can conclude 


-^7,.e 2 +rni a^)e-+b^e+ 

. V y 

(94) 

The highest eigenvalue of this matrix is the meta- 
canonical potential 


TOi -I- TO2 


E =-^7^cosh(/3ie-/2)-b ^ 
















where we have introduced 
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7 = - ^ 7;y sinh 



mi — TO 2 
2 


is 


(96) 


^ cosh[(aa - a^)e_] 

Jot — ! - 5 

y 7^ + J2uy 7 ^ 7 ^-e(“--“.')<^-/+(&.+fc. 0 E+ 


The meta-canonical potential T provides the statistics 
of ja the energy current flowing from the ath reservoir 
toward the system and of fa the activity induced by the 
ath mechanism. From direct derivation of F with respect 
to Oq, , ha or ruz the energy current coming from reservoir 
a > 1 is 


, _ E,.e- 7 a 7 i.e('’-+'’“)'^+sinh[(a„-a^)e_] 

Jot - j - 7 

y 7^ + Y.^,y 7,7,.eK-“.')^-+(&.+V)^+ 


the activity for the transitions induced by mechanism a 


and the occupation of state z is 


1 , (^ 2.1 - '^^. 2 ) 7/2 

Pz — a , — 


(99) 


where 7 of Eq. (96) is evaluated in m = c(a,2a), and 
taking & = 2 a to obtain the mean values of j, / and p 
in the NESS with affinity 2a. Deriving once more with 
respect to Ua', ha’ or mz' leads to the symmetric response 
matrix, see Eq. (78). The left and right eigenvectors of 
K associated to the eigenvalue F are respectively tt and 
r = • p. We find for the two-level model 


7^g(/3i/2-a,.)e_-sf,„£+ 7,,e/5i<=-/2 _ mi -f F ’ 

_-TOi -fF_ 

E,. 7!^e/5i<=-/2 _ rrii -f F ’ 

(Ejy (El, -I- Ei, 71 , 6 ^ 1 _|_ p) 

E„,„' 7„7i,-e(“>'-“-')E--s(b„-tV)£+ + -mi+vf 

(E„ - mi -F F) (Ei, + Y^ - mi -f F) 

Ei,,„' + [Yu 7i^e/5i^-/2 -mi+vf 


( 100 ) 

( 101 ) 

( 102 ) 

(103) 


We can now illustrate the consistency of the theory: from 
Eqs. (100-103) and the product tt • r, we recover the 
NESS probability of Eq. (99) obtained from derivation 
of the meta-canonical potential; Eqs. (100-101) allow us 
to compute the drift u to get the current and activity of 
Eqs. (97-98) from Eqs. (25-26) knowing the NESS prob¬ 
ability. 

We turn to the discussion of the properties of the two- 
level system with y = 2 heat reservoirs in the light of 
Fig. 3 obtained from our analytic results. For simplicity, 
we chose bi = 0. We set the energy scale and the time 
scale taking respectively /3i = 1 and 71 = 1. Fig. (3a) 
shows that the meta-canonical potential is a symmetric 
function of the affinity 02 and is strictly convex. From 
this symmetry, one should not conclude that the energy 
current j 2 is an anti-symmetric function of 02 - Indeed, 
the energy current comes from the derivative of the meta- 
canonical potential with respect to 02 evaluated in m = c 
that has no particular symmetry when changing the sign 
of 02 . 

The absolute value of the energy current \j 2 \ and the 
activity /2 always increases with the absolute value of the 
affinity I 02 I at given dynamical bias 621 see Fig. (3b) and 
Fig. (3e). A decrease of |j 2 | with increasing I 02 I would 


mean that the system has negative response for some 
affinities. Such a behavior is not expected for a simple 
two-level model. Another general trend is that \j 2 \ and 
/2 increases with 62 - Indeed, a higher dynamical bias 
increases the value of the transition rates corresponding 
to = 2, if one has e+ > 0, see Eq. (91). Then, a high 
dynamical bias accelerates the dynamics associated to 
reservoir v = 2, whereas a small one slows it down, let¬ 
ting the reference dynamics associated to reservoir 1 / = 1 
dominates in the transition rate matrix. Therefore, in 
the limit of low dynamical bias with respect to the affin¬ 
ity, the system approaches the EQ state at temperature 
/3i, with current j 2 and activity /2 decreasing to zero. 

To represent the NE micro-canonical potential L, one 
has to focus on the statistics of some specific variables 
by contraction: this step consists in evaluating the NE 
micro-canonical potential at the mean value of the disre¬ 
garded variables, for instance /i, pi and p 2 in the case 
of Fig. (3d). However, it is much more convenient to ob¬ 
tain L{j 2 , 72 ) directly from a parametric plot of {j 2 , f 2 , L) 
with ( 02 , ^ 2 ) being the parameters and taking hi = 0. 
In this way, we have obtained Fig. (3d) showing the NE 
micro-canonical potential as a convex function of (j 2 , / 2 )- 
This function is undefined in the regions corresponding 
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Figure 3. (a) meta-canonical potential for various ( 02 , 7112 ) with 6 = 0, (b) energy current, (c) energy, (e) activity and (f) 

entropy as a function of the affinity 02 and the dynamical bias 62 . (d) NE microcanonical potentials for the energy current j 2 
and the activity /2 after a contraction on /i, pi and P 2 - Other parameters are 61 = 0, P 2 = Pi — 2 a 2 , 72 = 0.5, ei = 1 and 
£2 = 0.5. For all hgures, /3i = 1 set the energy scale and 71 = 1 the time scale. The variables a and 6 are in unit of l//3i, the 
variables j and / are in unit of 'yi/Pi, and finally L, F, and m are in unit of 71 . 


to low activities in comparison to the energy current. 
The explanation is that a current can only flow if some 
minimal activity holds, i.e. if the system changes state 
regularly enough. 

Finally, the system energy e and Shannon entropy s 
are, in our framework, functions of the affinity and dy¬ 
namical bias. We see in Fig. (3c) and (3f) that these 
functions have a very similar shape in a large area corre¬ 
sponding to the EQ limit. The dimensionless free energy 
of the reference system at temperature /?i is (/? = /3ie — s 
and should reach its minimum value for low affinity I 02 I 
or low dynamical bias 62 . There, since /3i = 1, the system 
energy and entropy differ only in the value of the dimen¬ 
sionless free energy of the EQ reference system. On the 
contrary, at high affinity |a 2 |, most of the time the sys¬ 
tem is either in energy state ei = 1 for positive 02 , or 
£2 = 0.5 for negative 02 . The system is driven to a state 
where the entropy is lower than at EQ and the NE mean 
energy is moved away from the EQ mean value for the 
reference process. 


VI. CONCLUSION 

In this paper, we have established that the asymptotic 
probability of the energy currents, the activities and the 


occupations in a NE process proceeds from the long-time 
statistics of the same variables at EQ. We have connected 
the affinities of the NE process, the dynamical biases and 
the escape-rate changes to constraints imposed on the EQ 
reference process, respectively on the energy currents, on 
the activities and on the occupations of each state. This 
connection is the analog of the ensemble equivalence be¬ 
tween the canonical and micro-canonical ensembles of EQ 
statistical physics for which the temperature of the heat 
reservoir is associated to an energy constraint. We have 
argued that the mapping between EQ and NE fluctua¬ 
tions allows us to distinguish the reduced set of variables 
which play a key role in the description of NESSs. 

Beyond the understanding of the structure of NE 
statistical physics, phenomenological and/or operational 
methods must be developed to compute the NE poten¬ 
tials of real complex systems. In this regard, it was shown 
that efficient algorithms exist to compute the scaled cu- 
mulants of currents [ 66 ] or to simulate samples of rare 
trajectories [13]. A promising technique for macroscopic 
systems relies on the saddle point approximation of a 
path integral producing the cumulant generating function 
[67] . This calculation leads to a dynamical problem with 
a small number of degrees of freedom compared to the 
original problem. Solving this dynamical problem seems 
easier than finding the highest eigenvalue of a large tilted 
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operator. 
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Appendix A: LDF of empirical occupation and 
transition probability 


We derive in this appendix the LDF of transition prob¬ 
ability and occupation /(a),p) for the EQ process. By 
definition, the probability that oj[z\ = uj and p[z\ = p 
when the system trajectories are generated by the EQ 
reference process is 


Pt{i^,p) = '^'Pk[z]Sp^pi^]6^^^[z], (Al) 

[d 


or equivalently when introducing the empirical transition 
probabilities and occupations 


A[z] = dr ^vMr) 

Jo / 


,1, _ 

^xy 

Py 


E E ^x,z{T+dT)^y,z{T)^v,i/{T) In ^ 


TG[0,t] y,x^v,i' 


kxyPy ’ 


= ^ E 


UJ. 


Py[z] [Ky-'-^] + ^xy[z\ In 


xy 


kxyPy. 

(A5) 


Then, using the Kronecker symbols in Eq. (A2) one can 
take p[z\ = p and <jj[z\ = ijJ to move the action out of the 
sum and write 


Pt{uj,p) = e (A6) 


The remaining sum over all paths [z] is the probability 
that p[z\ and u}[z] take their typical value since the path 
probability is generated by uj/p. We expect this proba¬ 
bility to converge to 1 in the long-time limit. Then, we 
get the asymptotic probability of transitions and occu¬ 
pations of Eqs. (10-11) in the main text. 


Appendix B: Pre-exponential factor for the 
asymptotic probability of energy currents, activities 
and occupations 


where is the sum over all path [z]. We have intro¬ 
duced the action 


A[z] = In 


Vk[z] 


(A3) 


and the path probabilities with given initial state z(0) 


Vk[z] = exp 



P<..j/p[z] = exp 



Notice that the second line is identical to the first line 
where the empirical transition rate matrices ijA\z\Ip\z\ 
replace the real EQ rate matrices hy. From these path 
probabilities, the action becomes 



+ E In 

rG[0,i] 


i/(t) 

^z(T-\-dT)z{T) 


z{T+dr)z{T)xzp) 


Following, we provide an asymptotic approximation of 
the long-time probability of the energy currents, activi¬ 
ties and occupations when the final state is known. The 
variables are defined by 


ju[J\ 


fAz] 

Px[z] 


^ ^ ^ {^z(r+dT) ^z{t)]^i',iz{t)i 

rG[0,t] 

7 E dz(x^dx^z(x\S„ 


^2(r+dr)2(r) ? 


TG[ 0 ,i] 


drS.^ 


(Bl) 

(B2) 

(B3) 


The corresponding generating function with given final 
state X is by definition 


gx{a,b,m) 



(B4) 


It satisfies the linear equation dg/dt = k ■ g with k the 
tilted operator defined in the main text in Eqs. (35-36). 
We now look for a long-time asymptotic approximation 
of g 


gx{a,b,m) = ^ ^ 


(B5) 


^■^)xvPy> (B6) 


^ \ '' ^r{a,b,m)t 

-A-oo ^ 

y 

cz er(“’'’™)V,(a,6,m)(7rt./), (B7) 


>-oo 


(A4) 
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where is the initial state probability. We remind that 
TT and r are respectively the left and right eigenvectors of 
K for the highest eigenvalue F. Using an asymptotic ap¬ 
proximation to compute the inverse Laplace transforma¬ 
tion of gx{a,b,m), one recovers the exponent appearing 
in Eq. (12). 

Then, the pre-exponential factor in Eq. (B7) must 
be evaluated in o = a(j,/,_p), b = b{j,f,p) and m = 

We now assume that / and p are related 
to each other via the physical system subspace con¬ 
straint, see table II. From Eq. (49), we find 6, c) = 
7r2;(0,0, 0) = I for all x, where the second equality stands 
from the fact that the left eigenvector of a Markov matrix 
has all its components equal to one. Then, tt^ • = I 

by normalization of and the right eigenvector of k 
in the physical system subspace is the NESS probability 
r{a,b,c) = p* for the dynamics with energy current j 
and activity /. This leads to the asymptotic probability 
of energy currents, activities and occupations when the 
final state at time t is x 

Pt{jJ,p*,x) (B8) 


Appendix C: Asymptotic equivalence of the NE 
process and the canonical process 


The path probability of the canonical process, with 
generator /C, is defined by exponentially weighting the 
path probability of the EQ reference process: 

This tilting procedure is sometime referred to as canoni¬ 
cal conditioning. We show in this section that the above 
canonical process is asymptotically equivalent to the NE 
process defined in Sec. Ill A. Erom Sec. Ill C, it is also 
equivalent to the driven process. The connection be¬ 
tween the driven process and the canonical process was 
first obtained in Refs. [33, 34] and studied in depth in 
Refs. [29, 31]. 

From the definition of the CGF in Eq. (29), we have 


VkIv] - ((^ 2 ) 


Since Eq. (33) is satisfied for all y, we can write it for 
any state ?/(r) along the trajectory [y] 

r -t- (^^) 

and upon integration over the time r, one finds 

tv = tm'^ ■ p[y] + f dr(\(^) - Ay(.r)), (C4) 

Jo 

since A = A from Eq. (55). Finally, Eqs. (C2-C4) leads 


to the asymptotic equivalence of the path probabilities 

' ^ 1 i y{r+dr)y{r) > 


rG[0,t] 


:_).oo J-J- y(r+dT)y(r) 


tG [0,i] 


n{y]- 

t—¥oo 


(C5) 

(C6) 


Hence, we have proved that the canonical process cor¬ 
responds at long time to a NE process that can be re¬ 
alized experimentally changing the temperatures of the 
heat reservoirs and the dynamical biases. 


Appendix D: Fluctuation Theorem 

The FT is an essential property of the stochastic en¬ 
tropy production [36, 68]. According to this theorem, 
a stochastic positive entropy production is exponentially 
more likely than the opposite entropy production, i.e. an 
entropy destruction. On average this implies a positive 
entropy production in agreement with the second law. 
Therefore, the FT is a probabilistic statement of the sec¬ 
ond law, and as such it is a very fundamental property 
of NE phenomena. It was first derived with a long-time 
approximation, but since the mean entropy production 
always increases, a FT should hold at all time [69, 70]. 
Because the entropy production may be appropriately 
defined using different NE variables such as work, heat 
or particle currents depending on the experimental setup, 
the ET has many faces [56, 71-78]. Generally, the joint 
probability distribution of a set of time anti-symmetric 
variables summing to entropy production will satisfy a 
FT [79]. In our case, a linear combination of the cur¬ 
rents gives the entropy production rate 

<7 = 2a^ ■ j. (Dl) 

Accordingly, the LDF and scaled CGF for the NE process 
have a ET symmetry. This symmetry strongly relies on 
local detailed balance, in other word, on the symmetry of 
transition rates. We already used the local detailed bal¬ 
ance to show the equivalence of EQ and NE fluctuations. 
We show in this appendix that the fluctuation theorem 
(FT) is a consequence of the mapping between EQ and 
NE fluctuations associated to the symmetric nature of 
energy-currents fluctuations at EQ. Using Eq. (47), we 
find 

r(— 2a — a, &, 77!.) = F(—2a — a -I- a, & -I- 6, tti -1- c) 

= F(a -I- a, 6 -I- 777 -I- c6 -I- 6,77!. -I- c) 

= r(a, 5,77i), (D2) 

where we have used the fact that EQ fluctuations 
are symmetric in the reversal of affinities r(a, 6, tti) = 
r(— a,5,m). Similarly, from Eq. (48), it is straightfor¬ 
ward to see that 

W, f,P) - L{-j, f,p) = -2a^ ■ j = -cr, (D3) 
since we have L(j, f,p) = L{-j, f,p). 
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